Optical qudit-type entanglement creation at long distances by means of small 

cross-Kerr nonlinearities 
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Entanglement represents an important resource for quantum information processing, but its gen- 
eration itself requires physical resources that are limited. We propose a scheme for generating a 
wide class of entangled qudit-type states of optical field modes at sites separated by noisy medium 
when only weak optical nonlinearities are available at both sites. The protocol is also based on 
exploiting a weak probe field, transmitted between the sites and used for generation of quantum 
correlations between two spatially separated field modes. The idea of probabilistic entanglement 
enhancement by measurement is discussed, and corresponding scheme for measuring the probe field 
state with linear optics and photodetectors not resolving photon numbers is proposed. It is shown 
that the protocol is applicable in the case when decoherence, limited efficiency and dark counts of 
photodetectors, and uncertainty of nonlinear coupling constants are present. 

PACS numbers; 42.50.Ex, 03.67.Bg, 42.50.Dv, 42.65.Hw. 



I. INTRODUCTION 

One of the most intriguing features of quantum me- 
chanical description of physical objects consists in pres- 
ence of quantum superpositions and especially of entan- 
gled states [li]. At the very beginning of the quantum 
theory development such states were found to possess 
quite counterintuitive (from the classical point of view) 
properties |2|, i^] . During the 20th century the attitude of 
physicists to such quantum states evolved from perceiv- 
ing them as either evidences of theory incompleteness or 
interesting but quite useless features of the world ,1,- 5] to 
understanding the opportunities provided by such quan- 
tum objects for solving numerous tasks of information 
processing (6l-[llj|. It has been shown that quantum no- 
cloning theorem [T2| . [Ts'l provides unconditional security 
of quantum cryptography protocols [1-H, IT^ - [l6| . while 
using entangled states of quantum register can lead to es- 
sential speed-up of solving several classes of complicated 
computational tasks [MJ, [11, [H [11 . 

From this point of view entanglement represents an im- 
portant resource for different information tasks. On the 
other hand, entanglement itself requires some physical re- 
sources for its generation and, thus, can be considered as 
an intermediate step on the way from physical devices to 
accomplishing tasks of information processing. Similarly 
to many other situations, the resources available nowa- 
days for entanglement generation are limited, and one of 
the most important problems is to find methods for ob- 
taining results as good as possible using as few resources 
as possible. 

Among systems, being promising for efficient entan- 
glement generation, optical field modes take their place 
due to possibility of long-distance transmission with rel- 
atively low decoherence and quite simple manipulation 
of the states with linear optics and photodetectors. Con- 
sidering optical entanglement generation, one can divide 
necessary elements into two classes: quite simple ele- 
ments, available "freely" (linear optical devices, photode- 



tectors not resolving photon number, classical optical 
states), and resources — all other elements, being quite 
challenging for construction. It is impossible to gener- 
ate entanglement with the first group of elements only 
(see e.g. Refs. fl9l - [2lj |). and, therefore, some resources 
are necessarily required. The main resource, used for 
entanglement generation by optical methods, is nonlin- 
earity, either measurement-induced [2^, !23| or provided 
by interaction between the field modes in some medium. 
The corresponding physical resources are special detec- 
tors (e.g. resolving photon number (22l. [23j) and nonlin- 
ear media respectively. Entangled states can be created 
also by means of linear optics from nonclassical states, 
which in their turn require nonlinearity for their genera- 
tion. 

Many quantum information processing and communi- 
cation tasks (quantum cryptography, distributed quan- 
tum computation, teleportation of quantum states) re- 
quire distribution of entanglement between parties, sep- 
arated by noisy medium. For such tasks a high-quality 
quantum channel also becomes an important physical re- 
source. So, the two main resources, required for gen- 
erating entanglement between distant states with opti- 
cal methods, are nonlinear interaction and a quantum 
channel. Therefore, efficient entanglement generation 
corresponds in this case to creating strongly entangled 
states with weak nonlinear interactions and noisy quan- 
tum channels (only such systems are available nowadays) . 

In present paper we consider weak local cross-Kerr in- 
teraction as a resource for creating nonlocal nonclassical 
states of spatially separated optical field modes. Cer- 
tain pro gres s has already been achieved in this field of 
research |24| - |31| , especially in the case of creating qubit- 
like entangled states. However, it remained quite chal- 
lenging to create more general classes of entangled states 
with such limited resources. We propose a protocol for 
creating wide class of qudit-type states (including entan- 
gled states) with arbitrary dimensionality in continuous 
variable (optical) system using weak cross-Kerr nonlin- 



earity (as the main physical resources) , linear beamsplit- 
ters, detectors not resolving photon numbers, and sources 
of coherent states. We show that entanglement of the 
states, created with our protocol, can be higher than 
unity (which is the limit for qubit-type states, that can be 
created by previously proposed methods) and, thus, our 
protocol provides more effective use of limited physical 
resources. 

This paper is organized as follows. In the next section 
we discuss the main ideas of entanglement generation be- 
tween distant sites when local nonlinear interaction and 
non-ideal quantum channel are available. Then the main 
operations of the proposed protocol and corresponding 
state transformations are presented. Section |IV] is de- 
voted to discussion of the peculiarities of projecting the 
"raw" weakly entangled system state onto strongly en- 
tangled desired final state. It is shown that parameters 
of the protocol are determined in the unique way by the 
desired final state, and corresponding relations are found. 
In Section El we demonstrate several applications of the 
protocol to creation of nonclassical and entangled states 
of optical field modes. In the last section we prove ap- 
plicability of the protocol for entanglement generation 
under realistic conditions by taking into account deco- 
herence, limited efficiency and dark counts of photode- 
tectors, and uncertainty of nonlinear coupling constants. 

II. ENTANGLEMENT GENERATION WITH 
LOCAL CROSS-KERR NONLINEARITY 

Cross-Kerr interaction itself can be used for generat- 
ing entangled states starting from uncorrelated states of 
a pair of quantum objects. Suppose a field mode c in 
a coherent state \"f)c interacts with another system (an- 
other optical field mode a [27l. IsoMssj or an atomic system 
[2^ [281 113] )■ Due to the interaction the phase of the co- 
herent state amplitude 7 of the mode c is shifted by the 
value, proportional to the number of excitations n of the 
system a: 

\n)a\l)c^\n)ahe'nc, (1) 

where x describes effective strength of the interaction. 
If the initial state of the object a is a superposition of 
states with different excitation numbers (e.g. a coherent 
state I a) a in the case of field mode), the final state of the 
considered system will be an entangled state, composed 
by pairwise combinations of coherent states with differ- 
ent phases of mode c and number states of the system a 
(Fig. ma)). 

However, experimentall y o bserved nonlinear interac- 
tions are quite weak [35l - l38l |: the effective nonlinear- 
ity strengthes, predicted in the most promising 4- level 
atomic systems with electromagnetically induced trans- 
parency on the basis of theoretical calculations and ex- 
perimental data, have the order of x ~ I0~^ 10~^ [stI - 
|42 | . The magnitude of phase space displacement caused 
by such cross-Kerr interaction is proportional to \ j\Xi a-nd 
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FIG. 1: (a) Entanglement generation by nonlinear interac- 
tion, (b) Entanglement generation between sites, separated 
by noisy medium. Only weak entanglement can be gener- 
ated, (c) Entanglement generation between sites, separated 
by noisy medium, with probabilistic entanglement enhance- 
ment. Strongly entangled states are available in the case of 
successful measurement outcome. 



in general case the final state can be weakly entangled. 
One of the solutions of the problem consists in effective 
nonlinearity enhancement by using intense fields c: for 
I7I ^ 1 the displacement magnitude can be large enough 
(ItIx ^ 1) even for small x [24-26, 34]. 

Entanglement generation gets more complicated if we 
take into account not only limited available nonlinear- 
ity, but also noisy medium between the sites (Fig.[TJb)). 
Then decoherence strongly limits maximal possible am- 
plitudes of the transmitted field c (]7] ^ 1), while the 
amplitude of field a or the number of excitations in the 
atomic system is also limited by losses in local storage. 
Therefore, this simple scheme is not applicable for gen- 
erating strongly entangled states when the sites are sep- 
arated by noisy media. Hence, under realistic conditions 
a more sophisticated scheme, including some kind of en- 
tanglement enhancement, is required. 

One way of obtaining strongly entangled quantum 
state is to implement entanglement distillation [43-46]. 
This approach requires storage of quite a large number 
of initial weakly entangled states and can be quite chal- 
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lenging. Another solution of the problem can be based on 
probabilistic entanglement enhancement [13, [H, [13, [Ml ■ 
For this purpose one can design certain measurement, 
carried out at Bob's site, with successful outcome trans- 
forming initial weakly entangled state into a strongly en- 
tangled state. Direct measurement of the state of field 
mode c makes this mode inaccessible for any further use. 
Measurements, implemented with linear optics and pho- 
todetectors not resolving photon number on field modes, 
obtained by splitting the mode c, completely determine 
the state of the mode and destroy entanglement. There- 
fore, some kind of nonlinearity is required at Bob's site, 
too. It is quite natural to suppose that this nonlinear 
interaction is the same as the one at Alice's site (see e.g. 
Refs. [13,[li,[M[ll)- 

Several schemes, based on probabilistic entanglement 
enhancement by a measurement at Bob's site, have al- 
ready been proposed for entangling distantly separated 
field modes, when effective strength of nonlinear interac- 
tions X is equal to tt (strong nonlinear interaction) [27l |. 
and for entangling atomic qubits [1^ as well as for cre- 
ating qubit-type entangled states of optical field modes 
[30I Isij . when only small nonlinearity is available. 

In present paper we solve a more general task of de- 
signing the measurement scheme for creation of arbitrary 
qudit-type state (including entangled states) from a wide 
class of possible states in continuous variable system us- 
ing weak cross-Kerr nonlinearity, linear optical devices, 
detectors and sources of coherent states. The set of 
achievable final states of the field modes a and h has 
the form of a sum of phase-correlated pairs of coherent 
states of the modes: 

n 

where coefficients c„ are arbitrary and can be fixed in an 
appropriate way for obtaining the state, most useful for 
certain practical applications. 

III. OPERATIONS OF THE PROTOCOL 

We consider the following system (see Fig. Alice 
and Bob posses local field modes a and h referred to be- 
low as main field modes; the probe beam (ancillary mode) 
is denoted as mode c; mode d is a reference field, trans- 
mitted from Alice to Bob immediately before (or after) 
ancillary field for decreasing infiuence of dephasing in the 
quantum channel on the final states. 

Main field modes, the ancillary and the reference fields 
are prepared in coherent states (|a)a for Alice's mode a, 
|/3)b for Bob's mode h, \"f)c for the probe field c and \"f)d 
for the reference field d) . The initial state of the system is 
therefore uncorrelated one. Local cross-Kerr interaction 
of the modes a and c (with effective strength x — the 
phase of a coherent state of mode a increases by x radians 
per each photon in c) leads to generation of correlations 



between the number of photons in the mode c and the 
phase of the coherent state of the mode d. Then the 
field c is transmitted to Bob's site through the quantum 
channel. In this section we suppose for simplicity that 
all correlations are preserved by the channel. After local 
interaction of the modes h and c (the effective strength 
of the interaction is supposed in this section to be also 
equal to x): taking place afterwards, the phases of the 
main modes d and h become correlated with the number 
of photons in the mode c: 

l*l)a.c = E l«e^^")a \^)c , (3) 

n 

where Qn{l) = However, the correlations, 

generated in the system, are weak in a realistic case. For 
small ancillary field amplitudes I7P ^ 1, required for 
decreasing decoherence in the quantum channel i27l [28l . 
[3lj | , the entanglement between the states of the mode a 
possessed by Alice and the modes b and c possessed by 
Bob is much less than unity: 

E,c,h{xM'\7\')^l, (4) 

where h{x) — — a;log2 x — {1 — x) log2(l — x). 

As discussed in the previous section, the last stage 
of the protocol consists in detection of the probe beam 
state. This operation should lead (in the case of success- 
ful outcome) to transformation of the weakly correlated 
3-modes state into a strongly correlated 2-modes state of 
the form ©. 

The idea of obtaining the desired strongly correlated 
final state from the state Q is based on the following 
decomposition: 

where \'i'±)^i,^ denotes the part of the system state, or- 
thogonal to \y:>)ci and the vector \ip)c is uniquely defined 
as 

= const. ^«/Q;(7))|n>e- (6) 

n 

Thus, the transformation of the state ([3]) into the state 
^ for arbitrary coefficients c„ can be realized by success- 
ful discrimination of the state \ip)c from a set of ancillary 
mode states, containing \ip)c and a complete system of 
states, orthogonal to \ip)c- 

The most simple from theoretical point of view way 
of discriminating the state \ip)c is implementing projec- 
tive measurement described by operator = \ip)cc{^\ 
which satisfies the following relation: 

^J^l)a6c=l*/)afc®l^)- (7) 

However, implementation of such measurement with- 
out additional resources (e.g. nonlinear interaction) can 
be too complicated or impossible. Solution of this prob- 
lem is discussed in the next section, and the discrimi- 
nation technique, based on a special case of general class 
of POVM measurements — "elimination" measurements, 
— is proposed. 
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FIG. 2: Scheme of entanglement generation: Alice prepares entangled state of the field modes a and c; then she sends the 
ancillary mode c and the reference mode d through the quantum channel to Bob; having correlated the states of the modes h and 
c, Bob measures the final state of the ancillary field mode (with the help of the reference mode d) and, in the case of successful 
outcome, announces Alice under the classical channel that entanglement was generated (otherwise the set of operations is 
repeated). 



IV. DISCRIMINATION TECHNIQUE 

A. "Elimination" measurements. General 
consideration 

An important fact that can be used for designing re- 
quired discrimination protocol is that coherent states rep- 
resent a natural basis for linear optical devices (the set 
of coherent states is closed under linear optical trans- 
formations). From this point of view, it is convenient to 
base the discrimination scheme on comparing the state of 
the probe beam c with certain coherent states, obtained 
by splitting the reference mode d. However, coherent 
states are nonorthogonal. It is this problem that leads to 
impossibility of implementing the discussed above pro- 
jective measurement in a general case with linear op- 
tics. An important set of measurements, implementable 
with linear beamsplitters and photodetectors for coher- 
ent states, is the set of "elimination" measurements (see 
e.g. Refs. [28', "47] and Fig. below) — the measure- 
ments with successful outcomes manifesting that the field 
mode state is not certain coherent state. 

The concept of measurements with outcomes, mani- 
festing that the input state is not certain fixed state, 
was proposed in Refs. [isl - lsol l as a part of unambiguous 
discrimination of two non-orthogonal states. The con- 
cept of "elimination" measurements was generalized in 
Ref. [iyj for the case of discrimination between N sym- 
metric coherent states on a circle on phase plane with 
linear optical elements and photodetectors. In these pa- 
pers the system was supposed to be prepared in one of 
the states from a fixed finite set (e.g. symmetric coherent 
states on a circle: jae^^'''/^), fc = 0, N — 1). Discrim- 
ination of one of the states (e.g. |a)) is equivalent to 
elimination of all the remaining states from the fixed set 
(in the considered example \ae'^^^^/^), k = 1, ...,N ~ 1). 

Considering ancillary mode states in our protocol, one 
faces a more general situation. A complete set of inde- 
pendent states is infinite for a field mode, and the mea- 
sured state does not necessarily coincide with certain ba- 



sis vector, but can represent a superposition of basis vec- 
tors. Therefore, complete description of state transfor- 
mations, occurring when elimination measurements are 
carried out, requires a more rigorous operator definition 
of the concept of elimination measurements. 

We define measurement, eliminating state \ip), as any 
POVM-type measurement (described by POVM {Ai}, 
^JAi]+Ai = 1) with at least one of the outcomes (re- 
ferred to in this paper as successful) being characterized 
by operator, denoted in further consideration as A|^) , 
with the property 



=0. 



(8) 



Having obtained this outcome, one can with certainty 
conclude that the measured state was not the state \ip). 
Considering a set of measurements eliminating states 
we will require (simultaneously with Eq. (|S]) for all 
vectors {ipj)) commutativity of operators, characterizing 
successful outcomes of the measurements: 



A 



0. 



(9) 



This condition means that successful elimination of the 
state Itpi) must not destroy the result of elimination of 
the state \4>j) for any pair of the states {tpi) and lipj) from 
the considered set. 

The conditions ^ , (O lead to the following important 
implication, useful for designing the required discrimi- 
nation scheme. Let the set {ji/jo), be a (non- 
orthogonal) basis of finite-dimensional Hilbert space and 
operators >l|^o), A|i/>k) describe successful outcomes 
of corresponding elimination measurements. Any state 
from the considered state space 



4>k\iP- 



(10) 



subjected to successful elimination of the states 
{ji/'r), IV'if)}; is transformed into the state 
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.A 



0o^|^Ak)---^1>Ai)I^o) 



(11) 
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where contribution of \ipo) only did not vanish. The fi- 
nal state of the system differs from the state \ipo), but 
it is a completely defined state and can be transformed 
unitarily into jV'o) (in our protocol the ancillary mode 
is finally discarded and this transformation is not neces- 
sary). Therefore, one can state that successful elimina- 
tion of all the basis vectors except lipo) leads to successful 
discrimination of the state \4'o)- 

Generalization of the definition to the case, when sev- 
eral operators ^l^"*) : M = 1: 2, correspond to successful 

elimination of the state (IV'j))- = Oj is straight- 

forward: Eq. © retains its form and must be valid for all 
used instead of single operator A^^.y The result, 
described by Eq. (fTTI) . also remains valid. 



B. State discrimination with "elimination" 
measurements. From infinite-dimensional to 
finite-dimensional space 

The described above technique can be used exactly for 
discrimination of the state \(p)c of the ancillary mode 
c in a special case of finite-dimensional space of input 
states. Suppose that the final state of the whole sys- 
tem (modes a, b and c) can be decomposed using finite 
number of ancillary field mode states, and these states 
together with the state \(p)c span {K + l)-dimensional 
space. Then we can choose K such independent vectors 
IV'i)c, j — ^,---,K, orthogonal to \<f)c, that the system 
{|V'i)c, \i'K}c, \ f}c} is complete in corresponding sub- 
space. As shown above, discrimination of the state \ip)c 
corresponds to successful elimination of the vectors \ipj)c, 

Such finite-dimensional case can be realized, for in- 
stance, when the nonlinearity strength is equal to x = 
2tt/N, where N is an integer In these special case 
the state of the modes a, b and c after the nonlinear inter- 
actions can be represented as a sum of N terms consisting 
of coherent state of the form |7e^'^'^/'^)c of the ancillary 
mode c and corresponding entangled state of modes a and 
b. Discrimination of a fixed state \je^^''°^-'^)c (by elimi- 
nation of iV — 1 coherent states |7e^'^'^/'^)c, k ^ ko) maps 
initial weakly-entangled state onto strongly-entangled fi- 
nal state of the modes a and b. 

In general infinite-dimensional case, however, the state 
of the system cannot be decomposed in such a way and 
the described discrimination technique can be applied ap- 
proximately. For the considered state ([3|) and for small 
ancillary field amplitudes I7I <C 1 high accuracy of ex- 
ploiting the discrimination technique can be achieved by 
restricting consideration by a space of states with lim- 
ited photon numbers. Suppose the expression ([2]) for the 
desired final state has K + 1 nonzero terms, i.e. c„ — 
for n > K. Then Eq. © for the state \ip)c also has 
K + 1 nonzero terms and this state belongs to {K + 1)- 
dimensional subspace of states with limited photons num- 
bers n < K. According to Eq. ([3]), the probability of 



presence of more than K photons in the mode c is pro- 
portional to and is small for the considered sys- 
tem. Therefore, the state ([3|) with probability close to 
unity belongs to the same {K + l)-dimensional subspace 
as the state \(p): 



K 



(12) 



-\5^[''^)abc 
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where II I iJ^'^^"'^') abc 1 1 = O (|7|"'^+"'^) . If we construct a mea- 
surement, which leads to correct discrimination of the 
state \lp)c for the subspace, spanned by the state vectors 
with photon numbers n < K, the distance between the 
obtained state (for the ideal system) and the desired final 
state 1^/) will have the order not greater than 



W\ ^ 112 

abc 
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\2{K+1)) 
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(13) 



where pK is the probability of successful generation of the 
desired final state. Therefore, if the probability pK has 
the order at least O (|7p^), the error of exploiting the 
approximate discrimination technique will have the order 
O (I7P) which is much smaller than the decrease of the 
final state fidelity caused by non-ideality of the quantum 
channel. Calculations below show that this condition is 
fulfilled for the considered system. 

Thus, for generating the state |^'/)a6 composed by the 
sum oi K + 1 phase correlated pairs of coherent states 
of modes a and &, one needs to construct a scheme for 
elimination of K vectors |V'j)c, j — 1, 



C. Coherent states as the basis set for 
discrimination 

The next step of designing the scheme for discrimi- 
nation of the state \(p)c is construction of the set of in- 
dependent states {|V'j)c} in the form suitable for real- 
izing elimination measurements with linear optical ele- 
ments and photodetectors. Realization of such measure- 
ment for coherent states is known (see e.g. Ref. (47|). 
therefore, it is desirable that the basis states \ipj)c be 
coherent states \jj)c with amplitudes jj. In this case 
the orthogonality conditions have the form ^ {(p \ 7-,)^ = 

X) Cn7"/ (^Qn{'y)Vrilj = 0. Then the coherent states am- 
plitudes 7j can be found as K roots of the K-th order 
equation 



K 

n=0 



(14) 



and are uniquely defined for any given set {c„}. These 
statements are correct in the nondegenerate case. 



Presence of degenerate roots of Eq. leads to linear 
dependence of the set {|7j)}- In this case additional vec- 
tors must be added to the set to provide completeness. 
These vectors can be constructed in the following way. If 
the root 7™ has the multiplicity > !> then for x = 
holds 



On the other hand, 



= 0, 



< s < L 



1. 



(15) 



(16) 



= ^(^|(c+)^|7™),el^'"l^/^ 

P^ -(|16 p prove orthogonality of the vectors 

; = (c'^Y l7m)c to \ip)c for s < /,„ - 1. Instead 
copies of the vector |7„ 



we obtain Z,„ vectors 
1, that are (i) independent and (ii) 



Eqs. 

l7i^^) 

of Im 

)c; S = 0, Ir, 

orthogonal to \(p)c- 

Thus, in general case the set of vectors \jm)c — 
(c"*") |7m)c, where s — 0, — 1 and index m enu- 

merates distinct roots of Eq. ([T^ . is the required set of 
K independent vectors, orthogonal to |(/?)c. This set of 
states can be constructed for any desired final state of 
the form Q with finite number of terms and is uniquely 
defined for a fixed set of coefficients {c„}. Therefore, 
Eq. (I14|) provides the unique solution of the problem of 
generating any final state of the form ^ with minimal 
exploited resources. 
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D. Implementation: nondegenerate case 

In the nondegenerate case discrimination of the state 
\ip)c is based on quite a well known technique of elim- 
inating coherent states (Fig. [S]). For example, for 
eliminating a single coherent state |7i)c the measured 
state is displaced in phase space by the magnitude —71 
(the displacement is described by operator Dd^ji) — 
exp (— 71C"'" -|- 7*c)). The displacement operator trans- 
forms coherent state |7i)c into the vacuum state, and 
detection of photons in the field mode after the displace- 
ment (Fig. ^a,)) manifests that the measured state was 
not |7i)c- It should be noted that the considered detec- 
tors need not resolve photon numbers or be 100% effi- 
cient. 

Elimination of a set of coherent states {|7j)c} can be 
carried out in a similar way by splitting the field mode c 
into K modes (Fig. Elb)). Each of the obtained modes 
is used for eliminating one of the coherent states |7j)c- 
In this case obtaining photocounts (" clicks" ) from all the 
detectors corresponds to successful outcomes of elimina- 
tion of all the vectors {|7j)c,i = 1, ---jK} and, thus, to 
successful discrimination of \f)c- 

The displacement operators can be effectively real- 
ized by mixing the field mode c with additional refer- 
ence modes c?i, dfc in corresponding coherent states 



FIG. 3; (a) Elimination measurement for a single coherent 
state |7i)c. Dc{—^i) is the operator of coherent displace- 
ment with the magnitude —71. (b) Scheme for elimination 
of a set of coherent states \^j)c, j ~ 1, ...A', (c) Elimina- 
tion measurement for a single coherent state |7i)c with the 
coherent displacement operator Dc{—^i) being effectively re- 
alized by mixing the mode c with the reference mode di in 
coherent state |7i)di at linear beamsplitter BSi. (d) Scheme 
for elimination of coherent states |7j)c, j ~ 1, ...K, with lin- 
ear beamsplitters {BSi, BSk-i', BS[, -BS^-i) and 
photodetectors(Di,..., Dk); reference coherent states for im- 
plementing necessary displacements in phase space are ob- 
tained by splitting reference field d and applying phase shifts 

(j>l, cj)K-l- 



l7i)di) \lK)dK (prepared by splitting the reference 
mode d and applying additional phase shifts) at linear 
beamsplitters (Fig. [Sfc)). For example, for elimination 
of a single coherent state |7i)c one mixes the mode c 
with a single reference mode di at linear beamsplitter 
BSi with transmittance Ti = cos^ 6*1 (Fig. [31[c)). If the 
amplitude of the reference mode coherent state equals 
71 = —iji tan 9i and the measured state of the mode c is 
\jx)c, after mixing at the beamsplitter the modes will be 
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in the state | cos • (7^ + 71 tan^ 6*1 ) ) ^ | i sin 9i • (7^; - 71 )) rfi ■ 
The obtained state of the mode c can be used for some 
further operations, while the mode di appears just in the 
state, required for implementing elimination measure- 
ment: the amplitude of the initial coherent state \'yx)c 
is displaced by magnitude —71. Additional phase fac- 
tor i is irrelevant for detection of presence of photons 
in the mode, and decrease of the field intensity by the 
factor sin^ 6*1 = 1 — Ti influences only the probability of 
detecting photons and, thus, of obtaining successful elim- 
ination outcome. In an important case, when the mode c 
is discarded after the elimination measurement, the efH- 
ciency of the measurement can be improved by requiring 
Ti = 5 <^ \ (then sin^ 9i ~ 1 and success probability is 
approximately the same as for the scheme in Fig. [Sja)). 

The scheme suitable for implementing elimination 
measurements in general case of K coherent states 
{|7j)cjJ = ^t--tK} is shown in Fig. [3Ud). If the mea- 
sured state of the mode c is coherent state \'^x)c^ we re- 
quire the states of the reference modes dj after mixing 
with the mode c at beamsplitters to be coherent states 
1*9(72; ~ lj))d.j , where q is certain coefficient independent 
of the mode number j. These states correspond to dis- 
placement of the amplitude "fx of the measured coherent 
state by magnitudes and can be used for carrying out 
corresponding elimination measurements. 

Transmittances Tj of the beamsplitters BSj in this 
case are defined in the unique way by the above re- 
quirement of obtaining correct coherent displacements. 
As in the case of single coherent state |7i)c, we require 
the transmittance of the last beam splitter to be small: 
Tk = S <^ 1 (the mode c is discarded after implementing 
elimination measurements). The condition of dividing 
the amplitude of measured coherent state into equal 
parts q between the modes dj leads to the following sys- 
tem of equations for the beamsplitter transmittances (we 
define parameters dj by Tj = cos^ 9j ) : 



auitli — q, 

sin 02 cos6'i — g, 

sin 9k cos 9k- 1 ■ ■■■ ■ cos 9i = q. 



(17) 



Solving these equations together with the requirement 
Tk = S, one finds the following expressions for the trans- 
mittances of the beamsplitters: 

^ (K-j-l)(l-J) + l ^ K-j 

{K~j)il-S) + 1 ^K+l-j' '^'^> 

and coefficient q: 



q= [K 



1 - S 



-1/2 



1/Vk. 



(19) 



The amplitudes of the reference modes coherent states 
7m are also defined in the unique way by the require- 
ment of obtaining correct coherent displacements in im- 



plemented elimination measurements. The following re- 
current system of equation can be obtained: 



7j cos9j + i'jj-i sui9j — —iqjj, 
7j = Ij-i cos9j + i^j sm9j, 
70 = 0, 



(20) 



where jj is the amplitude of reference coherent state, 
mixed to the mode c by first j beamsplitters. The solu- 
tion of this system of equations is 



7i 



iq 



cos 9, 



{7j +sin2 0,(71 -f...+7j-i)} 



(21) 



As stated above, coherent states of reference mode dj 
with the required amplitudes 7, can be obtained by split- 
ting coherent state 7 of the mode d (Fig. [3]Jd)). The 
transmittances Tj = cos 9'j of the beamsplitters BS'j and 
the phase shifts are solutions of the following system 
of equations: 



7j = I cos 9i 



cos f 



smf 



I cos ( 



J = l,...,/^-l 

•COS Q'k-i1- 



e'^^7. 



(22) 



The system contains K complex equations for 2K real 
variables r{, T^_^, 0i, 4>k-i, R.e7, Im7. The 
solution of the equations is quite cumbersome, and we 
will provide it in explicit form only for certain special 
cases, discussed below. 



E. Implementation: degenerate case 

In the degenerate case we need to eliminate not only 
coherent states, but also the states of the form \^m)c ~ 
(£+)" |7m)c, created from coherent states by adding fixed 
number of photons (photons-added coherent states — 
PACS). We will show that elimination of PACS can be 
done exactly in the same way as elimination of coherent 
states |7j)c- 

For designing a scheme, suitable for elimination of the 
coherent state with one added photon, e.g. c'^|7i)c, it is 
useful to notice that when a beam with a single-photon 
excitation is split into two parts by a beamsplitter, the 
excitation can be detected in one of the two beams, 
but not in both of the beams simultaneously. If the 
mode c, the state of which is measured, is initially in 
the state c"^|7i)c, after splitting the mode into two parts 
(e.g. modes ci and 62) we obtain a superposition state 
7f(c^ + cJ)|7i/V2)ci|7i/n/2)c2, where the added pho- 
ton can be found in one of the modes, but never in the 
two modes simultaneously. Therefore, for all the terms 
of the superposition at least one of the modes is in co- 
herent state |7i/\/2). If we implement measurements, 
eliminating the state |7i/-\/2), for both of the modes, 
two successful outcomes can never be obtained, if the 
measured state of the mode c was c^|7i)c. Therefore, 
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two simultaneously obtained successful outcomes of the 
elimination measurements correspond to elimination of 
the PACS with one added photon c^\-fi)c- It should be 
noted that exactly the same scheme would be obtained, if 
we tried to eliminate two coherent states |7i)c and \"f2)c 
with equal amplitudes 71 — 72 by the method, suitable 
for nondegenerate case. 



Elimination of a PACS with s added photons (e.g. 

l7i'*^)c) can be carried out in a similar way, taking into ac- 
count that when a mode with s photons is split into s + l 
parts, it is impossible to detect a photon in each of the 
s + l modes simultaneously. Therefore, for elimination of 
the state |7j;'''*)c one can splits the mode c into s+l modes. 
When s photons, added to coherent state |7i)c accord- 
ing to the definition of the state |7j'''')c, are distributed 
between s + l modes, at least one mode appear in coher- 
ent state |7i/Vs + 1) without added photons. Therefore, 
s + l successful outcomes of elimination of coherent state 
\ji/Vs+ 1) for the s + l modes cannot be obtained si- 
multaneously if the initial state of the mode c is the state 

17^^-')^ (or any of the states |7i)c, h[^^)c, l7i''~^^)c 
with lesser numbers of added photons). Thus, obtaining 
successful outcomes of the s + l measurements, eliminat- 
ing coherent state |7i/Vs + 1), corresponds to elimina- 
tion of the state |7i'''')c, as well as of the states |7i)c, 
|7{;^'')c, 17^'' "'"■')c- As in the case, discussed in the pre- 
vious paragraph, exactly the same scheme could be used 
for elimination of the set of coherent states |7i)c, \ls)c 
with equal amplitudes 71 = 72 = ... = 7s if the case were 
considered as nondegenerate. 



The full algorithm of designing discrimination scheme 
for degenerate case can be summarized as follows. At 
first, K roots 7j of Eq. (fT4)) are found. Then the set 
of coherent states \jj)c is constructed. If the root 7™ 
has multiphcity Im > 1, "copies" of the coherent 
state |7m)c are replaced by Im independent states |7m)c, 
\jm^)c, |7m™ ^^)c- In the obtained set of indepen- 
dent states single coherent states are eliminated by the 
method, discussed in the previous subsection. Sets of the 

states \jm)ci \jm')c, IjrrT ^'')c with different numbers 
of photons, added to the same coherent state, are elimi- 
nated by carrying out measurements, eliminating coher- 
ent state \j„i) (with amplitude, decreased by splitting), 
for Im modes, obtained after splitting the mode c. There- 
fore, Im "copies" of the root 7^, appearing in the list of 
roots of Eq. (jl4p . correspond in the final discrimination 
scheme to ^m-fold elimination of the coherent state \ jrn)c- 
It means that one need not make any difference between 
degenerate and nondegenerate roots of Eq. ([Tl)) . eliminat- 
ing coherent states \"fj)c as many times, as they appear in 
the list of roots of Eq. (fTi]) . Thus, the scheme, designed 
for discrimination of the state \ip)c in nondegenerate case 
(Fig-ISd)), is also suitable for degenerate case. 



F. Implementation: mathematical description 

Quite interesting result of applicability of the same 
scheme for both nondegenerate and degenerate case can 
be given more rigorous mathematical proof on the ba- 
sis of operator definition of elimination measurements 
(Eqs. dED, ®). 

As shown in Appendix |^ (see Eq. (IA17P ). transforma- 
tion of the system density matrix in the case of successful 
outcomes of the K measurements, eliminating coherent 
states |7j)c, (photocounts obtained from all the detectors 
Dj) has the form 



[out] 
Pabc 



in-), 
he "< 



(23) 



where density matrices p'^^^ and p^™*^ describe the sys- 
tem state before and after implementing elimination mea- 
surements respectively; superoperator M (Eq. (IA9p ) de- 
scribes the part of system state transformation, which 
does not depend on measurement outcomes; the sets of 



operators {A. 



17.) c' "J 



— 1,2,...} correspond to successful 



elimination of coherent states 



i[,";/j7,)c = 0, n, = l,2,.. 



(24) 



(25) 



and satisfy Eq. 

The scheme is apparently suitable for nondegenerate 
case, and one needs to shows that it also can be used 
when some roots of Eq. (|14p are degenerate, i.e. that 
successful outcome of " Z^-fold elimination" of the state 
l7m)c, corresponding to the root 7^ with multiplicity l^m 
(or, in other words, successful elimination of Z,„ "copies" 
of the state |7m)c) corresponds to elimination of the 

states |7„i)c, \lm)c, |7m"~^^)c- 

The successful result of "Zm-fold elimination" of the 
state |7m)c is described by operator 



A 



("1) oc Cr- -7 Y 



(26) 



where r = 



Til 



+ 71;^ > Im- The expression 

{c-JmY {c+y^(c+ + -l,y{c~JmY 

contains powers of operator (c — jm) not less than r—s > 
— s > 1 for s — 0,...,Z„i — 1. Therefore, operator 
(1^51) corresponds to elimination of the states |7m''), s = 

0,-lm~l-- 



oc 



CX (C - 7m)'"'' |7m) = for S = 0, Im " 1, 



(27) 



which proves the conclusion made in the previous sub- 
section. 
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G. Final state for successful and "semi-successful" 
results of discrimination 

For obtaining the expression for the final state of 
the main field modes a and b after discrimination of 
the state \ip)c by elimination measurements, it is con- 
venient to introduce "phase-shifting" operator Fab — 
exp (^ix (^a^a + (it describes change of the main 

field modes state after cross-Kerr interaction with the 
mode c possessing 1 photon) and to represent state 
(described by Eq. ([3])) using operator of coherent dis- 
placement with operator-type argument: 

|^'l)a6c = (^Fabl) \a)a\l3)b\0)c, (28) 
with the following property: 

cDc (^Fabj) = FablDc [Pabl) , (29) 

which leads to significant simplification of Eq. (j23l) for 
the final state density matrix. As shown in Appendix \K\ 
(Eqs. (|Il9]) . (|X20)) '). the final state of main field modes 
a and b after elimination of all the states {|7j)c} and 
subsequent discarding the ancillary mode c is described 
by density matrix 

p..=9^"i7ri*'/L.,(*}i+o(i7r+^), (30) 

where 

1 ^ 1 

--T.^-^ab\-)aW)b = -\^f)ab- 
n—0 

The obtained expression means that the distance between 
the desired final state and the state pab, generated by the 
scheme, has the order O (I7P) and is small enough to be 
neglected when nonideality of the system is taken into 
account. 

It should be noted that in certain cases the final 
states, generated when "clicks" were obtained not from 
all the detectors, can also be useful ("semi-successful" re- 
sults). Such states are described by expressions, similar 
to Eq. (1311) but without multipliers, corresponding to the 
detectors (with numbers ni, n2, ■■■) that did not produce 
"clicks": 

m^ni ,nr>,... 

(32) 

This expression can be decomposed in the form, simi- 
lar to the desired final state ^ but with lower possible 
degree of entanglement: 

|vl/(ni, n2, ...))a6 = E ^nin,,n,, ...) |ae*>^")^ |/3e^^"), . 

n 

(33) 



E.g. for the case of absence of only one photocount the 
number of terms equals to K (instead oi K + 1) and the 
coefficients c„ can be found as 

n 

where n = 0, K — 1. 

In the next section we provide several examples of final 
states that can be generated by the protocol for successful 
and "semi-successful" discrimination outcomes. 



V. EXAMPLES 
A. Superpositions with correlated photon numbers 

As the first example of possible applications of the pro- 
tocol to nonclassical states generation we consider cre- 
ation of a superposition of states of modes a and b with 
correlated photon numbers. We show that such super- 
positions arise quite naturally in our protocol and then 
use them to illustrate general formalism, developed in 
Section HVl 

As discussed above (see Eq. ([1])), cross-Kerr interaction 
correlates photon number of one of the interacting modes 
with the phase of coherent state of the other mode. The 
state |^i)afcc, obtained after cross-Kerr interaction of the 
main modes a and b with the ancillary mode c, can be 
considered either as a superposition, where phases of co- 
herent states of the modes a and 6 are proportional to 
the number of photons in the mode c (Eq. ([3])), or alter- 
natively as a superposition, where the phase of coherent 
state of the mode c is proportional to the total number 
of photons in the modes a and b. The latter interpre- 
tation of the state \'^i)abc implies that discrimination of 
coherent state \'ye^^^)c of the ancillary mode c fixes the 
total number of photons in modes a and b to be equal to 
n. The final state |^'/)a6 in this case is a superposition 
of Fock states of the modes a and b with the number of 
photons in each mode varying form to n and the total 
number of photons being equal to n for each term. 

Mathematically this statement can be proved in the 
following way. The state \"^i)abc can be decomposed in 
the form: 

00 

\^i)abc = E mn)}ab\ie^nc, (35) 

ri=0 

where 

n 

{a)\m) a\n — m) b (36) 

is a superposition of states of the modes a and b with fixed 
total number of photons (equal to n); function Qmict) 
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is defined as Qjn{a) = -^e l"lV2^ (5™(a)Q„_m(a) = 

Qn{oi)^^ (™) ; assume for simplicity that a = /3. 

Suppose that only K terms are significant in the su- 
perposition (|55|): |Q„(a)| <C 1 for n > -ftT, 



I* 



K 

Ijabc ~ 

r!=0 



<J>(n))a6|7e''^' 



(37) 



Then successful outcome of elimination of K coherent 
states {|7)c, |7e'^^"-'^)c, |7e*^("+^))c, |7e*^^)c} 
with subsequent discarding of the ancillary mode c 
transforms the state \^i)abc into the following state of 
the modes a and b with correlated photon numbers, 
described above: 



|n — m)b X 
x{l + O(QK+i(a)/0„(a))}, 



(38) 



Generation of states of the form p8p can be described 
by general formalism, developed in Section IIVI For this 
purpose we find coefficients c„, for which the general fi- 
nal state |5'/)ab (Eq. (H])) is equivalent for the desired 
final state \^{n,K))ab (Eq. ([38])). Then coherent states 
amplitudes 7j can be found by solving Eq. (I14p . and ex- 
pressions for the parameters of discrimination scheme can 
be derived. 

The general expression Eq. ([2]) for the final state of the 
modes a and b can be transformed to the following form: 



coefficients c„ must be equal to cq = e'^, Ci = — 1 — e*-^, 
C2 = 1 (for unnormalized state). Amplitudes of the co- 
herent states |7j)c, eliminated in discrimination scheme, 
are equal to 71 = 7 and 72 = ^e^''^ in this case. Ac- 
cording to Eq. (dH]), transmittances of the beamsplitters 
55*1 and BS2 are equal approximately to Ti « 1/2 (for 
(5^1) and T2 = 5. The amplitudes of the reference 
coherent states, defined by Eq. (PT|) . are 71 « —i^i and 
72 ~ — * (71 +72) /(a/2(5). Solving the system of equa- 
tions (1121), one finds 0i = -x/2, T[ w 1-^72, 7 ~ -il2- 
For coefficients c„(ni) (Eq. ([M]) ). characterizing final 
state in the case of "semi-successful" outcomes of dis- 
crimination when the desired final state for successful 
outcome is described by Eq. (jH]), one obtains the fol- 
lowing expressions: co(l) = e*^, ci(l) = —1, co(2) = 1, 
Co (2) = —1. The final state, generated when "click" was 
obtained from detector D2, is approximately a vacuum 
state: 



\^[l))ab=\Q)a\0)b + O{\a\^), 



(43) 



while the state, generated when "click" was obtained 
from detector Di, belongs to the class of states, described 
by Eq. 

m2))ab = '"^"''^"t,'^^'-'"^^ + o{\a\) ^ m, i)u, 



V2 



and, therefore, can be useful for certain applications. 



(44) 



s=Q (n=Q 



by decomposing coherent states \ae^^'^)a and \f3e'^^")i, in 
terms of Fock states, where states \^{s))ab are defined by 
Eq. dMl). 

In order to obtain |5'/)a6 — \^{s, K))ab, the coefficients 
c„ must satisfy the following system of equations: 



K 



n=0 



c„e*^'''" = for s' = 0, 1, s - 1, s -f 1, K. (40) 



Before solving this system, it is useful to compare it with 
Eq. (fH|) for the amplitudes and to notice, that if co- 
efficients c„ satisfy Eq. (PP)) , K complex numbers je^^^ , 
s' = 0, 1, s— 1, s-|- 1, K , apparently represent the K 
roots of Eq. p4p . Then, coefficients c„ are defined in the 
unique way (except for overall normalization constant) 
by the complete system of roots and are equal to 



CO = ckHc''''' , CK-i^CKY.e'''''. (41) 

s' s' 

For example, if the desired final state is the following 



\^j)ab^m,2))ab^ 
|0),|2)fc + V2|l)a|l)6+|2)„|0)fc 
2V2 



+ 0(\a\ 



(42) 



B. Maximally entangled states for protocols with 
fixed number of detectors 

Another group of examples represents states of the 
form Eq. ([2]) with maximal entanglement, which is pos- 
sible for a scheme with fixed number of photodetectors 
K (and, therefore, with fixed number of terms in the 
expression ([2]) for the final state |4'/)ab). 

In the most simple case of schemes with 1 detector 
[K = 1) the coefficients c„ in Eq. ([2]), maximizing the 
final entanglement, can be found analytically as 



Co = 1, 



ci = — exp (~i (|c 



sm x) 



(45) 



where for simplicity we consider unnormalized final state. 
Additional condition, required for maximization of entan- 
glement in this case, is |a|^ = |/3|^ (for simplicity we will 
assume without loss of generality that a = /3) . 

The final state of the system is approximately a Bell 
state 

l*+)a5 = (l+)''l~)''+|-)'^l+)'')/^' (46) 

where {|-l-)a,b, |— )a.f)} is the orthonormal basis for the 
states of the modes a and b: |±)a,6 ~ {W)a,b ± 
g-z|a|^ smx|Q,g»x^^^^)^ Therefore, the protocol can be used 
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for generating qubit-type quantum states with maximal 
entanglement (equal to 1), possible for qubit systems. 
The set of states {|7m)c} consists of the only state 



|7i)c with the amplitude 71 = 76 



by Eq. ([Til) . The transmittance of the beam splitter BSi 
(Fig. ^c)) equals 6, and the amplitude of the reference 
coherent state is 71 — —171 (1 — 6) /6. 

For the scheme with two detectors {K — 2) the coeffi- 
cients c„, providing maximal entanglement, can be found 
analytically for systems with x ^ 1 and |a| = |/3| in two 
limiting cases: |apx^ 1 (low distinguishability of main 
field modes coherent states with and without phase shift 
equal to X- I {ae'^^ | a) | « 1) and |apx^ > 1 (high dis- 
tinguishability: I (ae*-^ I a) I ^ 1). 

For ^ 1 the final entanglement is maximal for 



ci 



-2 1 



Co = 1, 



(47) 



C2 = e-^^l"l 



The final state of the system is 

\ui)a\U'i)b + \/2|u2)a|'"2>6 + \u5)a\ui)b 



I* 



flab 



, (48) 



where 



\u2)a = g2(|a)ae'l"l''^ " \ae^'^) ae-'^^^^"^) , (49) 
\u^)a = <?3(|a)ae'l"l'^ + \ae^'^)ae-'\"\"^ - 

- go|ae'^)a), 

are orthonormal basis vectors for the mode a and the 



basis vectors 



for the mode b are defined in a simi- 



lar way (with a being replaced by coefficients qi are 
determined by the condition of orthonormality of the sys- 
tem of basis vectors. The final state possesses en- 
tanglement E — 3/2, which is higher than the maximal 
value, achievable for a pair of qubits. 

In this case the coherent states |7i)c and |72)c, ex- 
ploited in the detection scheme, possess close amplitudes 

2i|a| Transmittances 



71 



^2 = {l±iV2\ 



I«|V) 



76- 



of the beamsplitters -65*1 and BS2 are equal approxi- 
mately to Ti w 1/2 (for J < 1) and T2 ^ 5; ampli- 
tudes of the reference coherent states are 71 ss — *7i 
and 72 ~ — « (71 +72) /(a/S^)- The parameters of dis- 
crimination scheme, defined by Eq. (j22p . are equal to 
01 = \/2|a|x, r{ « 1 - ,572, 7 « -172. 

In the opposite limiting case |apx^ ^ 1 (such condi- 
tion is satisfied simultaneously with x ^ 1 fii^ the sys- 
tems with intense fields a and h) entanglement reaches 
the bound for 3-level system E = log2 3 ~ 1.58 when the 



parameters are 

Co = 1, 

The final state in this case has the form 



(50) 



\ui)a\u:i)b + \u2)a\u2)b + \uz)a\ui)b 

V3 



(51) 



where basis vectors \ui)a,b are defined by Eq. ((49)) above. 
Amplitudes of coherent states I71) and I72) are equal to 
71^2 = (l±i%/3)e2*l"l'x/2. The transmittances of the 
beamsplitters BSi and BS2 and the amplitudes 71, 72 of 
reference modes coherent states are the same functions of 
the amplitudes 71 and 72 as in the previously discussed 
limiting case. Solving the system of equations (0^ . one 
can find 0i = 7r/3, T{ w 1 - 5'^ /2, 7 « -^72. 

For intermediate values of distinguishability (|apx^ ^ 
1), as well as for greater numbers of detectors {K > 3), 
optimal coefficients c„ can be found numerically. The val- 
ues of maximally possible entanglement for schemes with 
fixed number of detectors are shown in Fig. (solid lines). 
One can see that the maximal possible value of entangle- 
ment grows with increase of the detectors number K , and 
in certain cases it may be considered as more effective use 
of fixed resources (nonlinear interaction, quantum chan- 
nel) than can be achieved in schemes with qubit-type 
entanglement. 

Fig. m also illustrates entanglement of final states, gen- 
erated by these schemes for "semi-successful" outcomes 
(obtaining photocounts from lesser number of detectors; 
dashed, dotted and dot-dashed lines in Fig. |4|). These 
final states possess non-zero entanglement and in certain 
cases can also be useful for solving information process- 
ing tasks. For example, the scheme with K — 3 detec- 
tors, optimized for obtaining maximal entanglement in 
the case of successful discrimination outcome, in the case 
of two absent photocounts can also produce maximally 
entangled state from the space of states |^'/)a6 with 2 
non-zero terms (line 3f in Fig. |4]). 



VI. NONIDEAL SYSTEM 

A. Considered types of nonideality 

In previous sections we assumed for simplicity that the 
quantum channel and photodetectors are ideal. In order 
to prove realizability of the proposed entanglement gener- 
ation method in realistic situations we discuss influence 
of system nonideality on the fidelity of obtaining final 
state. 

In real system decoherence and dephasing accompany 
all the stages of the protocol: implementation of nonlin- 
ear interaction, transmission of the probe beam through 
the quantum channel, storage of main modes a and b 
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FIG. 4: Entanglement of the final states, generated by 
the protocol, versus coherent states distinguishability |apx'^- 
Solid lines: maximal entanglement, possible for generation in 
the scheme with fixed number of photodetectors K: K = 1 
(1), K = 2 (2), K = 3 (3). Dashed lines (2a, 2b): entangle- 
ment of the final state, generated by the scheme with K — 2 
detectors with optimized coefficients c„ in "semi- successful" 
case (one photocount instead of two; entanglement does not 
depend on the number of the detector, which "clicked"). Dot- 
dashed lines (3a, 3b, 3c): final state entanglement for K — 3 
and 1 absent photocount. Dotted lines (3d, 3e, 3f; line 3f 
coincides with the solid line 1): final entanglement for K — 3 
and 2 absent photocounts. 



time, and its nonideality is supposed to effect the 
final state fidelity negligibly; dephasing in the res- 
onators is assumed to be small due to controllable 
laboratory conditions); 

(iv) decoherence and dephasing in the quantum chan- 
nel; 

(v) limited efficiency (probability of detecting a pho- 
ton, present in the field mode, equals A < 1) and 
dark counts of photodetectors (obtaining photo- 
count with probability ^ when the mode is in the 
vacuum state). 

For describing nonideality of the system we consider 
four stages of the protocol separately and find (super- 
)operators, describing difference between the states, ob- 
tained in ideal and nonideal systems. It should be noted 
that operator, transforming one fixed state into another 
fixed state (ideal state into nonideal one in the consid- 
ered case), is not defined in the unique way: its action on 
states from orthogonal space can be arbitrary. In further 
consideration we try to choose operators, acting on the 
main modes only (but not on the ancillary one), from 
sets of equivalent operators, transforming ideal density 
matrix into the nonideal one. Then the final state after 
implementing all the stages of the protocol is expected to 
be presented in the form of certain superoperator, acting 
on the ideal final state |vl'/)ab of the main modes a and 
b. 



in local resonators. All these factors can be taken into 
account by solving corresponding master equations "s^. 
However, for the system, considered in our paper, influ- 
ence of some of the factors on the fidelity of final state 
generation is supposed to be small. Therefore, for the 
purpose of simplifying further description, we take into 
account only the following factors that can limit applica- 
bility of our protocol: 

(i) decoherence of modes a and b during cross-Kerr in- 
teraction with probe beam c (decoherence of the 
mode c is assumed to have negligible effect due 
to much smaller amplitude of the probe beam: 
I7I ^ |a|,|/3|; qualitatively, decoherence of the 
mode c during cross-Kerr interaction infiuences the 
final state in the same way as decoherence of this 
mode during transmission through the quantum 
channel, but is weaker); 

(ii) inaccuracies in the nonlinearity values of the used 
Kerr media (we assume that effective strengthes of 
nonlinear interactions carried out by Alice and Bob 
are equal to Xac = X + ^Xac and Xbc ^ X + ^Xbc 
respectively and differ from the value x used in the 
scheme optimization); 

(iii) decoherence caused by nonideality of local res- 
onator at Alice's site (Bob's resonator is not used 
for storing part of an entangled state for a long 



B. Nonideal cross-Kerr interaction at Alice's site 

The first stage of the protocol is cross-Kerr interaction 
of modes a and c. We describe this nonlinear interaction 
in nonideal case by the following master equation: 



—p = ix [a'^ ac'^ c, p] + KiL{a)p, (52) 



where L{X)p 



Xp,X 



X,pX 



As discussed above 



(item (i) in the list of nonideality types), the only kind of 
nonideality taken into account by Eq. (j52p is decoherence 
of the mode d. 

Assuming that the duration of interaction is A<i, one 
can characterize the nonlinear interaction by effective 
strength Xac = x'^ti and relative loses rate Ai — 
g2rtiAti _ 1 (Ai (/g - /)//, where Iq and / are beam 
intensities before and after the interaction). 

In the ideal case {ki = 0) state transformation due to 
discussed cross-Kerr interaction is described by unitary 
operator 



Uac = exp(ixaca"''ac+c). 



(53) 



The initial uncorrelated state |a)a|/3)f,|7)c is transformed, 
therefore, into the superposition, where the phase of co- 
herent state of the mode a is proportional to the number 



13 



of photons in the mode c: 

oo 

?7ac|a)a|/3)6|7)c = 5]Q„(7)|ae 



6" 



(54) 



n=0 



Transformation of the system state in the nonideal case 
can be found by solving the master equation Eq. ([5^ 

with the initial condition p(0) — P^ahc^ where 

(55) 

is the initial uncorrelated density matrix (initial ampli- 
tudes a'''^ /3''*'\ 7'"' of coherent states in the nonideal 
case must be larger that the amplitudes a, /?, 7 that are 
expected to characterize final state). 

The solution can be found by representing the density 
matrix of the modes a and c in the form 



"1 ,"2 

X a{t)e'^''"' ) ( a(i)e*'^'*"^ 



\ni )Jn2\ 



(56) 



which is preserved during evolution. Substituting this 
decomposition into Eq. (j52p and solving the resulting sys- 
tem of differential equations, one can obtain the following 
expressions for the quantities a{Ati) and p„j„2(Ati) at 
the end of the considered stage of the protocol: 



p„,„,(Ati) « Q„,(7(°^)Q;,(7^°^)exp{|a(0)|2 x 
X Ai {ixac{ni -n2)~ xlci^i - "2)^) } , 



(57) 



(58) 



where smallness of nonideality is assumed for simplicity 
of derived expressions. The exponential factor in Eq. ([55)) 
describes influence of Kerr medium nonideality on the 
state, obtained after the interaction. For characterizing 
transition from the ideal state Eq. ([M)) to the nonideal 
one the following superoperator, acting in the state space 
of the mode a, can be chosen from the class of equivalent 
operators, describing this state transformation: 

Mi{a): |a(°)e*xacni ^^^^ „(o)g»xacn2 1 ^ 

^ |Q,'e»Xacni ^^^^Q,'g»xaon2| .exp||a(")p x (59) 

X Ai(^^ixac{ni - 112) - ^Xac("i - "2)^) }• 

This superoperator adds small phase shift to the coeffi- 
cients before coherent states of the mode a and decreases 
non-diagonal elements of the density matrix. 

Another factor, which influences fidelity of the final 
state generation but is not connected with nonideality 
of the Kerr medium itself, is deviation of the nonlinear- 
ity effective strength Xac from its expected value x, used 



during optimization of the discrimination scheme param- 
eters. This factor can be accounted for by introducing 
superoperator 



/ aa \ I 



(60) 



which provides additional phase shift to coherent states 
of the mode a. 



Storage of the mode a in nonideal local 
resonator at Alice's site 



The second stage of the protocol consists in transmis- 
sion of the ancillary field from Alice to Bob. At the 
same time the mode a, already correlated with the ancil- 
lary mode c, is stored at Alice's site. These two modes 
interact with the environment independently, and corre- 
sponding kinds of nonideality are considered separately. 

Decoherence of the mode a in nonideal local resonator 
is described by the following master equation: 



K2L{a)p. 



(61) 



If duration of this stage is equal to At2, relative losses 
rate equals to A2 = e^''^^*^ - 1. 

The solution of Eq. (j61l) can be found in the way, sim- 
ilar to the one used for the previous stage of the pro- 
tocol. The amplitude of coherent states of the mode d 
a{At2) = a'e~'^^*^ after this stage of the protocol must 
be equal to its final value a. 

The influence of decoherence in Alice's local resonator 
on the state of the system can be described by the super- 
operator 



M2{a): |aV^"<="i )aa( I 
^ [ae^xacni )^^(ae*^"<="^| -expllaf x 

X A2(iXac(ni - ns) - ^xld^i - ^2)^) |, 



(62) 



transforming the system state in the same way as Mi (a). 



D. Transmission of the ancillary field through 
nonideal quantum channel 

Interaction with the environment of the ancillary mode 
c during the second stage of the protocol is described by 
master equation 



— p = kL{c)p + rL(c+c)p, 



(63) 



where the first and the second terms describe decoher- 
ence and dephasing of the mode c respectively (see item 
(iv) in the list of nonideality factors). For characterizing 
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this type of system nonideality one can introduce rela- 



tive losses rate A 
A(j) ■ 



j2KAt2 _ I g^jjj mean phase error 



Due to commutativity of L(c) and L(c"'"c) (in the sense 
that L{c)L{c^c)p = L{c^c)L{c)p), the master equation 
Eq. (j63|) can be divided into two independent parts, de- 
scribing decoherence and dephasing. 

State transformation because of decoherence has the 
form: 



E 

ri=0 



(1-e 



— 2KAi2 Vi 



cV(c+)". 



(64) 



Dephasing of the mode c transforms the system state as 



,n2\ ^ \ni 



(65) 



As stated above, for certain simplification of further 
consideration it is useful to choose superoperators, act- 
ing on the main field modes only, from the set of equiv- 
alent superoperators, describing transition between ideal 
and nonideal case. For this purpose we take into ac- 
count that in the expression (|54p for the ideal system 
state number states \n)c of the ancillary mode appear in 
pairs with coherent states \ae'^^")a- Therefore, decrease 
of non-diagonal density matrix elements in the basis of 
Fock states of the mode c is equivalent to corresponding 
decrease of non-diagonal elements for the mode a in the 
basis of coherent states. In a similar way, discrete changes 
of photon number in the mode c, caused by energy losses 
in the quantum channel and described by Eq. (j64p . are 
equivalent to corresponding discrete changes of phase of 
the mode a. After certain mathematical calculations, one 
can show that the difference between ideal and nonideal 
states, caused by decoherence and dephasing of the mode 
c, can be described by superoperators 



and 



ji=0 



n 



Mf^(a): |ae'^'«^"i )aa("^*'' 



I aa \ 



iXo<:"2 I -A0-(ni-n2) 



(66) 



(67) 



respectively. 



E. Nonideal cross-Kerr interaction at Bob's site 

The third stage of the protocol is cross-Kerr interac- 
tion of modes b and c, described by the following master 
equation: 



b+h 



K3L{b)p, 



(68) 



and characterized by effective nonlinearity strength Xbc — 
x"Ai3 and relative losses rate A3 = e^'^^^*^ — 1 (we will 



assume for simplicity that A3 = Ai), where At3 is the 
duration of the interaction. 

In the ideal case this nonlinear interaction is described 
by operator 



Ubc = ex.p{ixbcb^bc^c), 



(69) 



with the obtained state UbcUac\ct)a\l3)b\j)c being equal 
to \^i)abc (see Eq. dS])). 

For nonideal system, master equation Eq. (I68p can be 
solved exactly in the same way as Eq. ([5^ . The superop- 
erator, describing transition between ideal and nonideal 
cases, has the form 



^ |/3e^X'>cni )^^(/?e*'^'"="^| •exp||/3(°)p x 
X Ai [^iXbc{ni - ^2) - ^Xbcini - "2)^) } 



(70) 



and, similarly to Mi{a), adds small phase shift to the 
coefficients before coherent states of the mode b and de- 
creases non-diagonal elements of the density matrix. 

Deviation of the nonlinearity eflfective strength Xbc 
from its expected value x can be accounted for by intro- 
ducing the following superoperator, describing additional 
phase shift to coherent states of the mode 6: 



Mx(6): |/3(°'e*x"i )fcfc( 
^ 1/3 



'bb^ 
^(0)g*X6cni 



(71) 



F. Nonideality of discrimination scheme due to 
limited efficiency and dark counts of photodetectors 

The last stage of the protocol is discrimination of the 
state \(p)c of the mode c at Bob's site. This stage in- 
cludes operations on the ancillary mode c only and does 
not influence directly modes a and b (their state is trans- 
formed due to previously generated correlations with the 
ancillary mode). All the superoperators Mi (a), Mp^(a), 

A'hid), M^^Ha), Mi^\a), Maib), M^{b), introduced for 
describing nonideality of the three preceding stages of 
the protocol, act on the state spaces of the mode a and 
b. Therefore, they must commute with any superoper- 
ators, characterizing the last stage of the protocol, and 
can be considered as acting after implementation of non- 
ideal discrimination measurement. For such considera- 
tion the input state of the discrimination scheme is the 
state |^'i)af,c, defined by Eq. 

For the ideal system, obtaining successful outcome of 
all the elimination measurements, followed by discard- 
ing of the mode c, transforms the input weakly entan- 
gled state |5'i)a6c into the desired final state \^f)ab (see 
Eqs. ([5n|) . (PTt ). Probability of this successful outcome 

equals p^'''"''^ = (^ItH^/IckP- 

Limited efficiency of detectors (described by the proba- 
bility A of registering photons — see item (v) in the list of 
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types of system nonideality) leads to decrease of probabil- 
ity of obtaining successful discrimination outcome, effec- 
tively reducing fraction q of the coherent state amplitude 
of the ancillary mode, interacting with photodetectors, 
by factor X: q Xq. Then one-run success probability 
for nonideal detection scheme is equal to 



Pk 



1 



\ck\ 



M 

K 



K 



1 



(72) 



Dark counts of photodetectors lead to mixing density 
matrices, characteristic to "semi-successful" outcomes, 
to the final density matrix, corresponding to successful 
elimination of all the states \^j)c- Then, according to 
Eqs. ((30)) ~(p4 | . successful discrimination of the state \ip)c 
by the scheme with nonideal photodetectors transforms 
the state |5'i)a6c into the following mixed state: 

P'ab=\^f),,J'^f\ + 



A|7P 
A2|7| 



\CK 



ni,n2 



*(ni,n2) 



ibab 



where j-th term corresponds to presence of j 
counts. 



(73) 



1 dark 



G. Final state in the nonideal case 

Summarizing the results concerning discussed types of 
system nonideality, we can express the density matrix of 
the final state of the modes a and 6 in the following form: 



( final) 
Pab 



Mi{b)M^{b)M^^\a)M. 



(74) 



This expression can be simplified by taking into ac- 
count that coherent states of the modes a and 6 posses 
correlated phases in all terms of the expression for the 
density matrix p'^^ (Eq. (|73l) ) and appear only in groups 
of the form c„ | ae*x") a \ ;3e*x" or 5„ (...) | ae^^" )^\(3e'^^}t. 
Therefore, superoperators Afi(d), Af2(a), Afj (a), Af3(6) 
act at the system state in the same way: they add small 
phase shifts to the coefficients c„ (or c„ (...)) and de- 
creases non-diagonal elements of the density matrix. Due 
to commutativity of the superoperators, they can be col- 
lected in a single superoperator 



Mo (a, 6) 



(75) 



transforming pairs of coefficients CnjC*^ of the state 
1^/ )( ■^^11 pairs of coefficients c„(ni, ^2, ...) 

with the similar meaning, defined by Eq. (j34p ) in the 
following way: 



„* pi';i(«l-n2)-')2("l-n2)^ 



(76) 



where 



ryi = -Ai {\aWac + m\bc) + |aPxacA2 (77) 



is a phase difference per photon and 



-\a\\lA2 (78) 



describes decay of non-diagonal elements of density ma- 
trix. 

The first term in the exponent of Eq. (|76)) corresponds 
to changing phase of coefficients c„ and can be compen- 
sated by corresponding changes in the detection scheme 
(by replacing 7^ by 7^6"*''^ in expressions for the scheme 
parameters). Therefore, only the second term of the ex- 
ponent is essential for estimation of the deviation of the 
nonideal final state from the ideal one. 

Finally, expression (|74p for the density matrix, ob- 
tained in for nonideal system, can be rewritten using the 
notations of Eq. ([75]) in the form 



^(/W) ^ M„(a,S)M(i'(a)Af^(a,S)pl„ 



(79) 



where M^{a,b) = M^{b)M^{a). 

In further consideration we assume that all types of 
nonideality, present in the system, are weak enough and 
the fidelity is close to unity (such systems are most use- 
ful from the practical point of view). Then the density 
matrix, defined by Eq. (j79p . is approximately equal to 



( final) 
Pab 



l*/)aba6(*/| + 

\^AMoia,b) + AA4'^\a) + 

+ AM,(a,S)}|*, + 



(80) 



A|7|^ 



where notation AMi = / — Mi is introduced for the su- 
peroperators; only the leading order of small parameters, 
characterizing nonideality of the system, is taken into ac- 
count. 



H. Fidelity of the final state generation 

According to the standard definition, the fidelity of 
generating the desired final state |^/)afc equals 



j^f \pi'r''\^f)ab 

ry, {final) 



(81) 



where is the actual final state of the main field 

modes a and S, defined by Eqs. (EH), ((50)) . 

Using Eqs. (|80| . (|8T|) and carrying out quite straight- 
forward mathematical calculations, one can obtain the 
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following expression for the fidelity of final state genera- 
tion: 



1 



J2 Tr{p_^AM,pif""') 



{AM.} 



} 



A|7l' 



|cK|^^Tr{Pj*K)U,(*K)|}, 



(82) 



where P_^-l - |*/ )^,^,( */|. 

For further consideration it is convenient 
to introduce unnormalized states |vE'j'^)ab = 

^^^g n'*c„|ae*'^")a|/3e*'<^")b, which are useful, for 
example, when exponent is decomposed in Eq. (|76p. 
Then it is quite easy to show that 



Tr 



(83) 

where we suppose that the term proportional to rji van- 
ishes due to correct phase compensation. 

Representing states of the form \ae^^'"'")a and 
|/3e^X''-")h as 



\ae ^ )n — e ^ \ae ^ ) 



(84) 



and decomposing exponents in power series in this ex- 
pression, one can simplify the term in the sum in Eq. (I82p , 
corresponding to deviations of nonlinearity strengthes 
from their expected value x'- 



Tr 



{^f I (Axaca+a + Axbcb+bf\ 



(85) 



Analytical expression for the remaining term in the 
sum can be found in two limiting cases: 



Tr 



(86) 



for la] X ^ 1 and 



Tr 



[p^Mi'\a)p^fr'^} ^ (87) 



for |apx^ > 1- 

Eqs. (1121), dMI), (ESI) -(113) provide the expression for the 
final state generation fidelity for any desired final state 
|*/)a6 described by Eq. 



I. Estimation of the system parameters, required 
for protocol implementation 

Further simplification of the derived above equations 
for the final state fidelity can be carried out in special 
cases, considered in SectionlVl We discuss the states with 



maximal entanglement, possible for a scheme with fixed 
number of detectors {K = 1 and K = 2). On the ba- 
sis of explicit expressions for the fidelity of such states 
generation we find conditions, that must be imposed on 
the system parameters in order to obtain sufficiently high 
fidelity of final state generation. 

For most of physical systems, suitable for implemen- 
tation of the protocol, assumptions that x ^ 1 and 
^ 1 are valid. In this paper we discuss analytical 
results obtained under these assumptions only (however, 
the opposite limiting case japx^ ^ 1 can also be de- 
scribed analytically). 

For maximally entangled state, generated by the 
scheme with K ^ 1 detector (Eqs. (|i5|). (HSl)). we ob- 
tain 



Tr 



and 



A|7l' 



Tr{p^AM,(a,S)pif"'^')} = 

^ 2\a?{eac + ebcf + {ea,-ebcf ^^^^ 
4 

"1 

c 



(90) 



2A|7P|a|V' 



where eac,bc — Axac,bc/x are relative inaccuracies of the 
nonlinear interaction strengthes; we assumed for simplic- 
ity that \a\ = and q = 1/Vk. 

For the scheme with K = 2 detectors and coefficients 
c„ described by Eq. (gZl) the form of Eqs. (I88l)-(l90l) re- 
mains the same, but expressions in the right hand side of 
Eqs. dm), dMl get numerical factor 2. 

The finally obtained expression for decrease of the final 
state fidelity contains 6 distinct terms, corresponding to 
different types of processes in the system: Eq. (1551) corre- 
sponds to 3 terms, proportional to A(/)^, Ai and A2 (see 
Eq. (UHl) for 772); Eq. dH]) and Eq. ^ provide expres- 
sions for the terms, describing nonlinearity strength inac- 
curacies and nonideality of photodetectors respectively; 
Eq. ([5S)) corresponds to decoherence of the mode c that 
can be described by effective discrete phase errors in the 
mode a. In order to estimate parameters values, suit- 
able for final state generation with sufficient fidelity, we 
require each of the discussed 6 terms to be not greater 
than some small value e, e <C 1 (then the fidelity will be 
not less than F > 1 — 6e). The obtained 6 inequalities can 
be used for finding 6 independent system parameters. 

For this purpose we divide parameters, describing the 
system, into four groups: 

(i) parameters A, \a\ characterize exploited "local" 
equipment (photodetectors, maximal field intensi- 
ties providing small decoherence during local op- 
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erations) and are supposed to be fixed by char- 
acteristics of existing equipment; in numerical es- 
timations we assume that dark count probability 
has the values C — 10~^ when detectors efficiency 
equals A = lO'^ and C = 10"^ for A = 10"! (such 
values can be achieved for InGaAs/InP photode- 
tectors fSU^); we also assume that laP ^ 10; 



(ii) parameters x, A^ac, ^Xhc Ai characterize local 
cross-Kerr interaction; we use the inequalities to 
estimate these parameters values and to find out 
whether such nonlinear interaction can be realized 
exp er inient ally ; 

(iii) parameters A, A0^, A2 characterize properties of 
the quantum channel and local resonator and de- 
termine maximal distance of entanglement genera- 
tion; 

(iv) the ancillary field amplitude I7I can be changed and 
is chosen so as to provide maximal success proba- 
bility for sufficiently high fidelity of the final state 
generation. 

Therefore, the discussed inequalities, providing suffi- 
cient fidelity, can be expressed in the following way (for 
K = \\ for K — 2 the conditions are the same except for 
numerical factor 1/2 in conditions 2, 3, 4, 6): 

A 2e2A 
A2 < 2e, 

Ai<-., 

k|2 < 



(91) 



|a|VA' 



2 a 



The first of the conditions limits acceptable losses in 
the quantum channel. This limitation is fundamental 
for optical methods of information processing when dark 
counts of photodetectors are present (see e.g. comments 
in Ref. [lal)- For instance, for the considered above pa- 
rameters and the desired fidelity F = 0.9 the maximal 
acceptable attenuation of the channel is limited by val- 
ues (20 28) dB, which for optical fiber with attenua- 
tion 0.20 dB/km correspond to maximal distances about 
Lmax ~ (100 ~ 140) km. However, it should be noted, 
that entanglement generation at such distances requires 
quite long storage of the field a in Alice's local resonator. 
The value of the resonator finesse, required for preserv- 
ing sufficient state fidelity and defined on the basis of 
the second condition of Eq. (PT|) . is about 10^^ in this 
case. Even greater values have already been predicted 
theoreti cally for crystalline whispering gallery mode res- 
onators [53|. Experimentally demonstrated high-quality 
resonators are characterized by values up to 10^ 10^^ 



[53l |5J| . Therefore, we believe that the protocol will be 
more suitable for efficient entanglement generation when 
the best available quantum channel connecting Alice's 
and Bob's sites is lossier than optical fiber. 

The third condition provides lower bound on the non- 
linearity value. For A0 ^ lO"'^ 10^^ and the consid- 
ered above parameters the minimal nonlinearity strength 
is Xmin lO^'^ 10"^. Such values has already been 
predicted in existing systems for the case of precise radi- 
ation focusing [37H43 |. It should be noted, that not only 
nonlinearity strength, but also acceptable signal attenu- 
ation during cross-Kerr interaction is limited (the fourth 
condition of Eq. (lOTl) leads to requirement Ai < 0.025). 

The fifth condition of Eq. (|9T|) limits maximal probe 
beam intensities and determines the maximal possible 
one-run success probability for the protocol, described 
by Eq. (1721) . For the protocol with K = 1 detector the 
success probability is sufficiently large for all losses values 
not exceeding the limit determined by the first condition 
of Eq. (Fig.jS]). For K — 2 generation of the desired 
final state can be implemented without too large number 
of ancillary field transmissions for losses not more than 
approximately 14 dB (which correspond to the distances 
up to 70 km in optical fiber). 

The maximal acceptable relative inaccuracies of non- 
linearity strengthes Eqc, £fcc, defined by the last condition 
of Eq. jnH), are equal to 0.09 for if = 1 and to 0.06 for 
K ^2. 

Such parameters values can be achieved in real sys- 
tems, and, therefore, our calculations prove applicabil- 
ity of the protocol for entanglement generation between 
sites, separated by lossy media, using contemporary ex- 
perimental equipment. 



VII. CONCLUSIONS 

To summarize, in the present work we have proposed 
a protocol for creating a wide class of qudit-type states 
(including entangled states) with arbitrary dimensional- 
ity in continuous variable system using weak cross-Kerr 
nonlinearity, linear beamsplitters, detectors not resolving 
photon numbers, and sources of coherent states. 

The method of entanglement generation is based on us- 
ing an ancillary field mode, transmitted from Alice's site 
to Bob's one through lossy quantum channel. Weak non- 
linear interaction of the mode with the main field modes 
possessed by Alice and Bob leads to creation of a weakly 
entangled 3-modes state. The main problem, solved in 
our work is designing a scheme for the ancillary mode 
state measuring leading to probabilistic entanglement en- 
hancement and transforming the "raw" weakly correlated 
state into highly entangled final one. The found POVM 
measurement is shown to be implementable with linear 
optics and photodetectors, not resolving photon num- 
bers, on the basis of elimination measurements. The 
equation, defining parameters of the detection scheme 
for a given desired final state in a unique way, is also 
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FIG. 5: (a) Dependence of one-run success probability pK on 
the relative losses rate A for the protocol with K photodetec- 
tors for fixed fidelity of final state generation F — 0.9: K = 1 

(1) , K — 2 (2). Dotted line corresponds to the limit of entan- 
glement generation because of dark counts of detectors, (b) 
Dependence of one-run success probability pk on the desired 
fidelity F for the protocol with K photodetectors and differ- 
ent losses rates A: = 1, A = 14 dB (1), A' = 1, A = 28 dB 

(2) , if = 2, A = 14 dB (3), A' = 2, A = 28 dB (4). 



derived in our paper. 

Our calculations prove applicability of the proposed 
protocol in the case of realistic photodetectors with lim- 
ited efficiency and dark counts, nonlinear Kerr interac- 
tion with decaying modes, uncertainty of coupling con- 
stants, and lossy quantum channel. It has been shown 



that the protocol can be used for creating quantum states 
with entanglement higher than unity and, therefore, in 
certain cases corresponds to more effective use of re- 
sources than can be achieved for protocols based on en- 
tangling qubit systems. The fidelity of final state gener- 
ation F — 0.9 can be achieved when a quantum channel 
with losses rate up to (20 -;- 28) dB is available (it corre- 
sponds to distances up to 140 km for optical fiber). Re- 
quired cross-Kerr nonlinearity is x > Xmin ^ 10^^-^ 10^^ 
and can be created using contemporary equipment. 
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Appendix A: Elimination of coherent states as 
POVM measurements 



Here we provide mathematical description of system 
state transformation when measurements, eliminating 
coherent states {|7j)c}, are implemented. The consid- 
ered measurement scheme is shown in Fig. EKd). The 
mode c is mixed at linear beamsplitters with K reference 
modes dj, prepared in coherent states \ jj)dj by splitting 

the initial coherent state of the mode d. Then photode- 
tectors Dj determine presence of photons in the modes 
dj. 

Let the state of modes a, b and c before elimination 
measurements be described by density matrix p^^^^ . Tak- 
ing into account that coherent states represent a natu- 
ral basis for describing linear optical transformations, it 
is convenient to decompose the input density matrix in 
terms of coherent states of the mode c 



(Al) 



where Pati'^x) is an operator- valued function of variable 
(acting as an operator on the modes a and 6), analo- 
gous to Glauber function of one-mode field. 

The state of the expanded system, composed by the 
main modes a and &, the ancillary mode c and refer- 
ence modes dj, after mixing the ancillary mode c with 
reference modes dj at beamsplitters is described by the 
following density matrix: 



J 



Pahcd,...dK = TcdK-Tcd^ {^Pabc ® )rfidi(7l| » - » \lK )dKd,AlK\j f+^- 



.T 



cdK ' 



(A2) 



where f^d^ = exp ^iOj (c+dj + cdj^ | are the unitary operators of field transformation by beamsplitters (the 
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quantity 9j is related to transmittance as Tj = cos^ 6j). 
Using Eq. (jAll) one can represent the density matrix 
in the form 

Pahcdi...dK — 

= /^a.(7.)l^(7.)),,,...,^,,,...,^(*(7.)M'7., ^^^^ 

where \^{-ix))cd^...dK = T,,dK■■■Tcd^\lx)c\ll)d^■■■\ll)dK■ 
T}le values of transmittances Tj of the beamsplitters 
BSj (Eq. (ITSl) ) and amplitudes 7^ of the reference modes 
dj (Eq. (j21[) ) are chosen in such a way, that amplitude of 
coherent state \'^x)c of the ancillary mode is split in equal 
parts between the modes dj, and each of the reference 
modes effectively undergoes coherent displacement — 7^ : 

\^{lx))cd^...dK = Wo{lx))c® ^^^^ 

® Vq ■ iix - 7i))rfi ■•■ \iq ■ {ix - iK))dK , 
I 



where 7g(7a;) is a linear function of the amplitude ^x of 
the ancillary mode coherent state, equal to 



7o(7.) = - Kq^ 7, + y qY^lK- (A5) 

Measuring presence of photons in the modes dj by the 
photodetectors can be described by a set of pairs of pro- 
jective operators 

^j- - 1 ) d, ( 1 and Pj+ = l-Pj-, ( A6) 

describing absence and presence of photons in corre- 
sponding field mode respectively. 

The final state of modes a, b and c after carrying out 
the measurements and discarding reference modes dj is 
described by density matrix 



{out) 
Pabc 



K 

PathxMhx) ),,(7o(7x)l ■ IlTrd, {Pj±Uq-hx-jj) >d^.d^.(*'7- (7. -7,)l} • A., (A7) 



i=i 



where the type of used projector ("-I-" or "-") depends on the obtained measurement outcome. This expression can 
be simplified using the following relation: 

Trd,{Pj±\iq-ilx-lj))a^aiiq-hx-lj)\] ^ ("j ' (7x " 7j)) = 



E 



(A8) 



where fl- = {0} and fi+ ^ {n \ n > 1} are the sets of photon numbers, corresponding to the measurement outcomes 
"-" (absence of photons detected) and "-f" (photocount obtained) respectively. 

Introducing superoperator, which describes measurement-invariant part of the ancillary mode state transformation 
by the definition 



M: |7. )Jlx\ ^ |7o(7x) >,,(7o(7x)l exp^-g^ ^ I7. - 7jf ) , 
one can transform Eq. (jA7p for the density matrix /o[j™*^ to the form 

Pa,{lx)\lx)Jlx\-\{ 



{out) 
Pabc 



Then one can define operators, which act as follows 



= —T^ - Ij) ' ■ \lx)c ^ (rix - 7j) ' \lx) 



and transform the expression (jA10[) for the final state density matrix in the following way: 



(A9) 



(AlO) 



(All) 



{ou€) 
Pabc 



M J j Pabilx) ■ Y BKn^...BinAlx ) Jlx\B+,,-B + ^, ■ d^Jx 



(A12) 
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Eq. (jA12p describes transformation of the system den- 
sity matrix by elimination measurements in the form, 
similar to the one, corresponding to POVM measure- 
ments. It should be noted, however, that, the standard 
normalization condition is satisfied only for complete 
state transformation by the measuring setup (including 
action of superoperator M and K operators 
Bkuk) rather than for single state elimination. 

Operators Bjnj for rij > correspond to the definition 
Eq. ([5]) of operators, describing successful elimination of 
the state \jj)c- 



Bjn,\l3)c 



ilx ~ IjT' \1t)c = 0, rij = 1,2,, 



(A13) 

Therefor, returning to the notations of Eqs. 

we can define operators AQnj{\lj) c), eliminating coher- 



ent state |7j)c, as 



(c-7j)"\ = 1,2,... (A14) 



These operators satisfy the conditions, provided by 
Eqs. dll), (HD: 



^l7,>c'^|7,>c 



= for all 



(A15) 



(A16) 



According to Eq. (|A12p . the final state of the system 
in the case of successful outcome of elimination of all the 
states {|7j)c} is described by density matrix 



If the input state is |4'i)a6c = Dc (^Fab^^ \a)a\f3)b\0)c (see Eqs. (pS)) ). the expression for the final state can be 



rewritten as 



TiiL-.TiK'- \ 7 / V 7/ uucuuL. y ^„ J ^ 

(A18) 



where Eqs. (|A14p . were taken into account. For 

small ancillary field amplitudes I7I ^ 1 the main contri- 
bution to the final state density matrix is made by the 
term with ni = ... = uk = 1 (the most probable case of 
successful elimination of the states {|7j)} corresponds to 



detection of exactly 1 photon by each of the detectors). 
After discarding mode c, the final state of modes a and 
b is described by density matrix 



Pab = Tr,pi;f = q'^'hn^'f )^,J^'f\ + O (|7r+^) , (A19) 

where 

\^)ab = (Pab - ^) ... (Pab - ^) \a)Mb = — E CnP:b\c^)a\l3)b = — |*/)a6; (A20) 

I 



we have taken into account that amplitudes {7^} are scribed by Eq. 
roots of Eq. (jl4l) : \^f)ab is the desired final state, de- 



[1] E. Schrodinger, Proc. Cambridge Philos. Soc. 31, 555 [2] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
(1935). 



21 



777 (1935). 

E. Sdirodinger, Naturwissenschaften 23, 807 (1935). [29 
N. Bohr, Phys. Rev. 48, 696 (1935). 

N. Bohr, in Albert Einstein: Philosopher-Scientist, Ed. [30 
P. A. Schilpp, Evanston, 1949, p. 200. 

C. H. Bennett and G. Brassard, in Proceedings of IEEE [31 
International Conference on Computers and Systems and 
Signal Processing, 175 , (Bangalore, India), (1984). [32 
A.K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 

C. H. Bennett, G. Brassard, and N.D. Mermin, Phys. [33 
Rev. Lett. 68, 557 (1992). 

D. Deutsch, in Proceedings of the Royal Society of Lon- [34 
don, Vol. A 425, 73, (1989). 

P.W. Shor, in Proc. of the 35th FOCS, 124, (1994). [35 
L. Grover, in Proc. the 28th ACM STOC, 212, (1996). 
W.K. Wootters and W.H. Zurek, Nature 299, 802 (1982). [36 

D. Dieks, Phys. Lett. A 92, 271 (1982). 
Ed. M. Pavicic, Quantum computation and quantum, [37 
communication: theory and experiments, Spriringer, [38 
2006. 

Eds. N. Cerf, G. Leuchs, and E. S. Polzik, Quantum In- [39 
formation with Continuous Variables of Atoms and Light, [40 
Imperial College Press, 2007. 

M. Dusek, N. Lutkenhaus, and M. Hendrych, in Progress [41 
m optics, Ed. E. Wolf, Elsevier, Vol. 49, 2006, p. 381. 
Eds. G. Leuchs and T. Beth, Quantum Information Pro- [42 
cessmg, WILEY- VCH GmbH & Co. KGaA, Weinheim, 
2003. [43 
S.Ya. Kilin, in Progress in optics, Ed. E. Wolf, Elsevier, 
Vol. 42, 2001, p. 3. [44 
M.S. Kim, W. Son, V. Buzek, and P. L. Knight, Phys. 
Rev. A 65, 032323 (2002). [45 
W. Xiang-Bin, Phys. Rev. A 66, 024303 (2002). 
J. K. Asboth, J. Calsamiglia, and H. Ritsch, Phys. Rev. [46 
Lett. 94, 173602 (2005). 
T.C. Ralph, A. Gilchrist, G.J. Milburn, W. J.Munro, [47 
et al., Phys. Rev. A 68, 042319 (2003). [ 

E. Knill, L. Laflamme, and G. J. Milburn, Nature 409, [49 
46 (2001). [50 
P. van Loock, T. D. Ladd, K. Sanaka, F. Yamaguchi, [51 
et al., Phys. Rev. Lett. 96, 240501 (2006). 

S.G.R. Louis, K. Nemoto, W.J. Munro, and T. P. [52 
Spiller, New Journal of Physics 9, 193 (2007). 
S.G.R. Louis, K. Nemoto, W.J. Munro, and T. P. [53 
Spiller, Phys. Rev. A 75, 042323 (2007). 
S.Ya. Kilin and A.B. Mikhalychev, in Proc. SPIE, [54 
Vol. 6726, 6726 3D, International Society for Optical En- 
gineering, (USA), (2007). 

P. vanLoock, N. Lutkenhaus, W. J. Munro. and K. 



Nemoto, Phys. Rev. A 78, 062319 (2008). 

P.P. Rohde, W.J. Munro, T.C. Ralph, P. vanLoock, 

et al., QIC 8, 0053 (2008). 

S.Ya. Kilin and A.B. Mikhalychev, Nonlinear Phenom- 
ena in Complex Systems 12, 150 (2009). 
S.Ya. Kilin and A.B. Mikhalychev, Optics and Spec- 
troscopy 108, 178 (2010). 

T. Tyc and N. Korolkova, New J. Phys. 10, 023041 
(2008). 

D. Mogilevtsev, T. Tyc, and N. Korolkova, Phys. Rev. 
A 79, 053832 (2009). 

T.D. Ladd, P. van Loock, K. Nemoto, W.J. Munro, 

et al.. New Journal of Physics 8, 184 (2006). 

J. Toulouse, Journal of Lightwave Technology 23, 3625 

(2005). 

J.H. Lee, T. Nagashima, T. Hasegawa, S. Ohara, et al.. 
Journal of Lightwave Technology 24, 22 (2006). 
H. Kang and Y. Zhu, Phys. Rev. Lett. 91, 093601 (2003). 
H. Kang, G. Hernandez, J. Zhang, and Y. Zhu, LANL 
Archive: quant-ph/0510036 (2005). 

H. Schmidt and A. Imamoglu, Opt. Lett. 21, 1936 (1996). 
M. D. Lukin and A. Imamoglu, Phys. Rev. Lett. 84, 1419 
(2000). 

G. F. Sinclair and N. Korolkova, Phys. Rev. A 76, 033803 
(2007). 

S. Li, X. Yang, X. Cao, C. Zhang, et al., LANL Archive: 
arXiv:0802.2136vl [quant-ph] (2008). 
C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, 
et al., Phys. Rev. Lett. 76, 722 (1996). 

C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. 
Wootters, Phys. Rev. A 54, 3824 (1996). 

M. Horodecki and P. Horodecki, Phys. Rev. A 59, 4206 
(1999). 

G. Alber, A. Delgado, N. Gisin, and I. Jex, J. of Phys. 
A 34, 8821 (2001). 

S.J. van Enk, Phys. Rev. A 66, 042313 (2002). 

I. D. Ivanovic, Phys. Lett. A. 123, 257 (1987). 

D. Dieks, Phys. Lett. A 126, 303 (1988). 
A. Peres, Phys. Lett. A 128, 19 (1988). 

D. Stucki, C. Barreiro, S. Fasel, J.-D. Gautier, et al., 
LANL Archive: arXiv:0809.5264vl [quant-ph] (2008). 
J. Zhang, P. Eraerds, N. Walenta, C. Barreiro, et al., 
LANL Archive: arXiv:1002.3240vl [quant-ph] (2010). 
I. S. Grudinin, V.S. Ilchenko, and L. Maleki, Phys. Rev. 
A 74, 063806 (2006). 

S.M. Spillane, T.J. Kippenberg, K.J. Vahala, K. W. 
Goh, et al., Phys. Rev. A 71, 013817 (2005). 



